# NOT RUN {
data(Pb)
# }
# NOT RUN {
#Default method:
set.seed(1023)
boot.Pb.default <- bootMSD(Pb$value, Pb$u) # Uses individual standard uncertainties
summary(boot.Pb.default)
#Method for MSD object:
msd.Pb<-msd(Pb$value, Pb$u) # Uses individual standard uncertainties
boot.Pb <- bootMSD(msd.Pb, B=5000)
#Increased replication compared to default
summary(boot.Pb)
# NOTE: The default summary gives individual observation p-values.
# To correct for multiple comparisons, apply
# a suitable p-value adjustment:
summary(boot.Pb, p.adjust="holm")
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab